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We present the thermal infrared (5-35 pm) spectrum of 956 Elisa as measured by the Spitzer Infrared 
Spectrograph ("IRS'’; Houck, J.R. et ai. [2004]. Astrophys. j. Suppi. 154, 18-24] together with new ground- 
based lightcurve data and near-iR spectra. From the visible lightcurve photometry, we determine a rota- 
tion period of 16.494 ± 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the 
visible absolute magnitude (Hv) at that rotational phase to be 1 2.58 ± 0.04. From radiometric analysis of 
the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 
5,3 ± 0.4 km with a visible albedo p v = 0.142 ± 0.022, significantly lower than that of the prototype V-type 
asteroid, 4 Vesta. (This corresponds to a radius of 5.2 + 0.4 km at lightcurve mean.) Analysis with the stan- 
dard thermal model (STM] results in a sub-solar temperature of 2923 ± 2.8 K and beaming parameter 
1.16 + 0.05, Thermophysicai modeling places a lower limit of 20 J nr* K“‘ s~i on the thermal inertia 
of the asteroid’s surface layer (if the surface is very smooth) but more likely values fall between 30 and 
ISO j m~ 2 K' s s~t depending on the sense of rotation. 

The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled 
thermal continuum, exhibits mineraiogically interpretable spectral features within the 9-12 pm restst- 
rahlen band, the 15-16.5 pm Si-O-Si stretching region, and the 16-25 pm reststrahlen region that are 
consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fail within the narrow 
compositional range Wo 2 -iEn 74 ± 2 Fs 2 4 - T . Spectral deconvolution of the 9-12 pm reststrahlen features 
indicates that up to «s20% olivine may also be present, suggesting an olivine-diogenite-iike mineralogy. 
The mid-lR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface 
of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the 
diogenite component. 

Analysis of new near-lR spectra of 956 Elisa with the Modified Gaussian Model (MOW; Sunshine, j.M., 
Pieters, C.M., Pratt, S,F. [1990], j, Geophys. Res, 95 (May), 6955-6966) results in two pyroxene composi- 
tions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene, High-Ca pyroxene is not evident in 
the mid-lR data, but may belong to a component that is underrepresented in the mid-IR spectrum either 
because of its spatial distribution on the asteroid or because of its particle size, High-Ca pyroxenes that 
occur as exsolution lamellae may also be more evident spectrally in the NfR than in the mid-IR, In any 
case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of dioge- 
nite-like composition, which has very rarely been observed among asteroids. 

Published by Elsevier Inc. 


1. Introduction 

The dynamical family of small basaltic asteroids associated with 
4 Vesta (“Vestoids”; Binzel and Xu, 1993) is considered to be the 
proximate source of the HED (howardite, eucrite, and diogenite) 
meteorites. The majority of these meteorites are thought to have 
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their ultimate origin on Vesta itself (McCord et al„ 1970: Consol- 
rrtagno and Drake, 1977). Hubble Space Telescope imaging of Vesta 
(Thomas et al„ 1997a) has revealed an impact basin large enough 
(D as 460 km) to be the source of all known Vestoids. Spectral stud- 
ies of Vestoids in the visible and near-lR have characterized their 
olivine and pyroxene compositions and revealed significant differ- 
ences between HED and Vestoid populations. For example, about 
25% of HED meteorites are diogenites. although to date, observa- 
tions of diogenite-like near-IR spectra among Vestoids have been 
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Fig. t. Proper elements of 956 Elisa in the context of 4 Vesta and its dynamical 
family. The small dots represent ail known asteroids; dark filled circles are Vesta’s 
dynamical family; and the large filled black dot represents 4 Vesta: The small filled 
triangles ate all spectroscopically confirmed V-types. 

very rare. A number of vestoids (Burbine et ai„ 2001 ; Duffard et al„ 
2004) have been characterized as having olivine-rich compositions 
long unrepresented among meteorites. (Recently, however, the 
meteorite MIL 03443 has been suggested (Mittlefehldt, 2008} to 
be a dunite sample from the HED parent body.) Similar N1R fea- 
tures suggest regional olivine exposures on Vesta itself (Gaffey, 
1997; Binzei et al„ 1997). 

The inner-main-belt Asteroid 956 Elisa was identified as a V- 
type asteroid by Florczak et al. (2002), who considered it to be out- 
side the limits of the Vesta dynamical family although it is rela- 
tively dose to the Vestoids in proper element space (Fig. 1). 
Subsequently, however, numerous other V-type asteroids were 
discovered in the inner-main-belt outside the boundaries of the 
Vesta family. A dynamical mechanism was proposed by which 
956 Elisa could have been transported to its present orbit from 
within the Vesta family (Carruba et al., 2005, 2007), and 956 Elisa’s 
orbit is now generally considered consistent with an origin on 4 
Vesta (Nesvomy et al., 2008). 

We present the thermal infrared (5-35 pm) spectrum of 956 
Elisa as measured by the Spitzer Infrared Spectrograph ("IRS": 
Houclr et al., 2004) together with new measurements of the visible 
lightcurve and near-infrared (0.8-2.5 pm) spectrum. The thermal 
flux spectrum, with constraints provided by the absolutely cali- 
brated visible lightcurve. provides new estimates of the radius, 
thermal inertia, and albedo. The lightcurve data also permit accu- 
rate determination of a rotation period for 956 Elisa, Mid-IR spec- 
tral emissivity features are also observed in the spectrum of Elisa: 
comparison with laboratory spectra of relevant samples, including 
HED meteorites, establishes that these features are diagnostic of 
mineralogy and permits the determination of a pyroxene composi- 
tion for 956 Elisa, The near-IR spectrum is used to independently 
estimate the surface composition of 956 Elisa, and the results of 
the two techniques are compared. 

2. Circumstances of the Spitzer observations 

Spectra of 956 Elisa were acquired on 2008-August-02 in the 
four low-resolution orders of the Spitzer IRS (Houck et aL, 2004). 
All observations were completed between UT 18:03:50 and UT 
18:26:26. During the observation, 956 Elisa was at R = 1.832 AH. 
A = 1.1 10 AH, and a phase angle of 29.32°. ! The acquisition mode 


1 Following convention, we use R to represent the asteroid-Sun distance and 4 to 
represent the asteroid-observer distance, (n this case, the observer is at the location of 
the Earth-trailing Spitzer observatory rather than on Earth. 
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Fig. 2. Spitzer IRS spectrum of 956 Elisa with standard thermal model (‘’STM") fit 
overplotted. The wavetength ranges of the four Spitzer low-resolution orders (S12, 
SL1 . LL2. IL1 ) are marked in blue . Data from the four orders were spliced together to 
assemble this composite spectrum. The scaling factors were: +1.12% (SL2/SL1), 
-0.84% (OZ/SLl), and -4.47% (LL1/SL1). Of these, the EX2/SL1 (14.15 pm) is the 
most uncertain due to an IRS instrumental artifact (Section 3.2) that sporadically 
produces excess flux in the 13.5-14.15 pm region of Sil . immediately shortward of 
the LL2 boundary. Thus, we have incorporated an additional 2.5% uncertainty in this 
value into our radius and temperature solutions. (For interpretation of the 
references to colour in this figure legend, the reader is referred to the web version 
of this article,) 

of the IRS was very similar to standard ground-based spectroscopy 
in that the asteroid was observed at two nod positions approxi- 
mately 1 /3 and 2/3 of the length of the slit. 

Much of the data reduction was carried out at the Spitzer Sci- 
ence Center (SSC) by pipeline processing. Data were returned to 
the observers as fiat-fielded, flux-calibrated spectral images. We 
removed background emission by subtracting frames from the 
two nod positions, corrected for rogue pixels using software sup- 
plied by the SSC (http://ssc.spitzer.caitech.edu/postbcd/irs- 
clean.html), and extracted the 1-D spectra using SPICE, another 
piece of software supplied by the SSC (http://sscspitzer.cal- 
tech.edu/postbcd/spice.html). For Elisa, we had obtained multiple 
cycles for each of the four modules at both nod positions. We aver- 
aged all cycles and both nod positions to produce a final spectrum 
for each module. The four modules were scaled relative to each 
other using flux densities in wavelength regions of overlap. The en- 
tire spectrum was then rescaled to the average flux of the four seg- 
ments, conserving the measured photometric flux. The scaling 
factors in ail cases fell within the 5% per order photometric flux 
uncertainties of the IRS. 

3. Thermal model fitting: radius and temperature 

in the 5-38 pm region, the spectrum of 956 Elisa is dominated 
by thermal emission. To model the thermal continuum of 956 Elisa 
(Fig. 2), we apply a variant of the Standard Thermal Model (“STM”: 
Lebofsky et al., 1986; Spencer etaL, 1989) to the IRS emission spec- 
tra of the asteroid. The STM begins with the assumption that the 
asteroid is a spherical greybody in equilibrium with the incoming 
solar radiation at every point on its surface: 

T(0)^T mBX coskO) (1) 

where 6 is the angular distance between any point on the sunlit side 
of the asteroid and the sub-solar point. 

The thermal flux from the surface is then calculated by integrat- 
ing the Planck function over the projected surface area of the ob- 
ject. For a spherical body, this can be expressed as 
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F(X) = ~ c(x) f B(X. T(0 )) cosf 0) sin(ff)dfl (2) 

A Jo 

where r is the asteroid radius, e(x) is the emissivity, and B(/.,T) is the 
Planck function. The two free parameters in this model, radius and 
maximum (sub-solar point) temperature, are varied to fit the mea- 
sured flux. In this relatively simple model, other factors such as al- 
bedo and beaming parameter are calculated directly from the radius 
and sub-solar point temperature. The particular strength of the IRS 
data for thermal modeling is the broad and continuous wavelength 
coverage from 5,2 to 38 pm that includes the thermal emission 
peak for asteroids. The temperature is tightly constrained by the 
shape of the spectrum and the size by the absolute flux level. 

3.J. Radius estimate 

While the procedure for estimating size is relatively straightfor- 
ward, several sources of uncertainty related to the absolute flux of 
the spectrum and some model parameters must be considered to 
assess the accuracy of the result. 

The formal error in the radius resulting from the fitting proce- 
dure based on the relative uncertainties of the Spitzer data points 
is quite small, ±0.1 km. Other sources of error dominate the overall 
uncertainty in our radius estimate. The most important are the 
absolute flux uncertainties in the Spitzer IRS data (±5% per order, 
thus ±2.5% for the average of the four orders), uncertainties in 
the scaling between the orders (+2.5%), and the thermal 1R phase 
correction. 

The thermal 1R phase correction can be made either with a 
phase coefficient which is constant with wavelength, established 
empirically as 0,01 ± 0.005 mag./' (Spencer et ai., 1989; Matson, 
1972); or by the NEATM method (Harris, 1998) in which the por- 
tion of the asteroid out of the field of view is eliminated directly 
from the STM flux calculation. In other words, the original STM 
integrates thermal flux at zero phase angle over the sunlit hemi- 
sphere and then decreases the flux by 0.01 mag./° to the phase an- 
gle of the observation. The NEATM, on the other hand, integrates 
the flux directly at the phase angle of the observation over the vis- 
ible hemisphere (T - 0 K on the nightside) as seen by the observer. 
The difference in final flux between the two phase curves depends 
strongly on phase angle and for moderate phase angles {<70°} 
peaks at 1 3-20% around 25-40°. 

Examination of the infrared phase curve of Ceres, for instance 
(Spencer, 1990) reveals that the NEATM provides a better match 
to the pre-opposition phase curve, but the 0.01 mag,/° is a better 
model for the post-opposition phase curve. Thus, neither model 
can be established to be universally more accurate. 

Because of the high phase angle (29.319°) of the Spitzer obser- 
vation, the ±0.005 mag./° cited uncertainty in the IR phase curve 
produces the largest component in the uncertainty of our radius 
calculation; +7.0%[-6.S% in radius. (At this angle, the NEATM pro- 
duces a phase correction equivalent to the lower range of the STM, 
0,005 mag./ 0 , so the STM/NEATM variation is included within this 
range.) 

As shown in Eq, (2), the integrated flux depends equally on r 2 
and bolometric emissivity. Uncertainties in the assumed emissivity 
will therefore propagate to the derived radius. The effective emis- 
stvity of a silicate asteroid with the surface temperature distribu- 
tion of 956 Elisa can be calculated based on the emissivity data 
in the ASU library (Christensen et ai., 2000): 

fi Fsm(x)€(X)dA 

€ effective — 1 77\T^ Wy 

h F S TM(A}dA 

Assuming that the asteroid's composition is 50-100% pyroxene and 
0-50% plagiociase, the effective emissivity will be 92.6% ± 1,7%. The 


corresponding contribution to the uncertainty in the radius esti- 
mate is 0.9%. 

Overall, we estimate that the projected radius of 956 Elisa at the 
time of the Spitzer observation was 5.3 + 0.4 km. 

3.2. Temperature estimate 

In the STM and NEATM, the temperature distribution of the 
asteroid's surface is a function of the maximum temperature (Eq. 
(1)). The maximum temperature from the thermal fit, T^*, is sub- 
ject to fewer sources of uncertainty than the radius estimate. 
Essentially, factors that change the shape of the measured or mod- 
eled spectrum will affect the temperature determination. The rela- 
tive (point-to-point) spectral uncertainty contributes only a 
negligible amount to the temperature uncertainty. An inherent 
wavelength dependence of the NEATM infrared phase curve is 
potentially responsible for a 1.1 K change in the temperature esti- 
mate; the uncertainties associated with splicing the four Spitzer 
IRS orders together may account for another 2.5 K. (Caution is 
especially required in considering the SL1/LL2 scaling, as the 
14.15 pm region is close to an artifact known as the “14-micron 
teardrop" [“IRS: Image and Spectral Features", http://ssc.spit- 
zer.caltech.edu/irs/features.html] that appears sporadically in SL1.) 

Including all of these uncertainties, we estimate the maximum 
temperature of 956 Elisa at the time of observation to have been 
T mix - 282.3 ± 2.8 K. 

The derived maximum temperature and size depend on the po- 
sition of Elisa, which is well known for all of our observations, but 
also contains information about the albedo, effective emissivity, 
thermal inertia, and surface roughness of the asteroid: 

_ 171 ~ A)(insolation)V 

T “ = w j (4; 

where A is the Bond albedo, a is the Stefan-Boltzmann constant, e is 
the effective emissivity, and ?/ is the "beaming parameter”, which 
represents a combination of the thermophysical properties of the 
surface (e.g. thermal inertia) and the surface roughness. 

In order to isolate these factors from one another, additional 
information such as reflected-light (visible or NtR) photometry is 
needed. The combination of visible absolute magnitude and mid- 
IR derived radius permit a solution for the albedo, thus enabling 
the beaming parameter and thermophysical properties of the 
asteroid to be further constrained. 

3.3. lightcurve and shape 

Visible wavelength lightcurves of 956 Elisa were measured to 
constrain its absolute magnitude, its rotational phase at the time 
of the Spitzer observations, and a first order estimation of its shape. 
These measurements were obtained at two observatories: the Uni- 
versity of Hawaii 2.2 m telescope and a 40 cm telescope at the 
Observatoire de la Cote d’Azur. The data from the Observatoire 
de la Cote d'Azur were generously provided by Matthieu Conjat 
and were obtained over the course of nine nights between 2008- 
June-28 and 2008-July-12. All of these images were obtained using 
a Johnson R-band filter. These data were presented and made avail- 
able through the CdR online data base of asteroid lightcurves 
(http://obswww.unige.ch/~behrend/page_cou.html). 

The UH 2.2 m observations were performed over the course of 
three nights from 2008-Juiy-26 to 2008-juiy-29. Images were ob- 
tained using an SDSS-r filter. Reduction of these data followed 
standard 1RAF procedures. The observing conditions on 2008- 
july-26 were photometric, which resulted in a mean error of 
0.014 on the measured magnitudes for that night. Though the mag- 
nitudes measured on the non-photometric nights were calibrated. 
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Fig. 3. Measured photometric lightcurve (in apparent SDSS-r band magnitude) of 
956 Eiisa. Data points from 12 separate nights of observations have been phase* 
folded onto a single period. The best fit lightcurve (solid line) displays a period of 
16,494 ±0.001 h and a peak-to-peak amplitude of 0.348 magnitude, The Spitzer 
observations viewed the asteroid at one of two possible rotational phases (shaded 
grey regions). These two possibilities (labeled PG for prograde rotation, and RG for 
retrograde rotation) exist due to the uncertainty in the orientation of Elisa's rotation 
axis. 


we did allow a small additive offset to be applied to these data 
when determining the best fit iigbtcurve. Fig. 3 shows the phase 
folded lightcurve for 956 Elisa with data from all 12 nights of 
observations. 

The measured magnitudes were fitted with a lightcurve of the 
form: 


r(t) = ]T a, sin [ 2 Jti + b, cos j ' 2rd — , (5) 

where r(r) is the r-band magnitude as a function of time. The coef- 
ficients ( di.bt ), phases , <^ 2i ) and period (P) were treated as free 

parameters in the fitting algorithm. As mentioned above, small 
additive offsets were applied to the non-photometric data. Individ- 
ual offsets were chosen for the set of data points from each night. 
These offsets were treated as free parameters and were restricted 
to within the photometric error bars. The best fit lightcurve was 
determined based on a minimization of the reduced chi-squared 
statistic. The fitted rotation period of the asteroid, P, was deter- 
mined to be 16.494 ± 0.001 h. The high precision in this measure- 
ment was facilitated by an observational baseline of more than a 
month and resulted in rotational phase errors at the time of the 
Spitzer observations that were insignificant (<1°). 

The orbital phase angle of Elisa (defined as the Sun-Elisa-obser- 
ver angle) as seen by the Earth-trailing Spitzer observatory was 
different from the orbital phase as viewed from Earth. The 
Earth-based light curves were measured at an orbital phase of 
~7“ whereas Spitzer viewed Elisa at a phase angle of 29.3°. 
This difference results in a degeneracy in the assumed rotational 
phase observed by Spitzer. This degeneracy is caused by the 
unconstrained orientation of Elisa's spin axis, if Elisa rotates in a 
prograde sense then Spitzer saw the surface of Elisa associated 
with a rotational phase of -340°. If Elisa's rotation is retrograde, 
Spitzer saw the surface of Elisa associated with a rotational phase 
of ~25° (shown as -385° in Fig. 3). This degeneracy can not be 
broken without knowledge of the orientation of Elisa's spin axis. 
Fortunately, the difference in magnitudes corresponding to these 
two rotational phases is small (Fig. 3). 

In the prograde case, Spitzer observed Elisa at an apparent mag- 
nitude of r- 13.75; in the retrograde case, Spitzer observed at an 
apparent magnitude of r= 13.73. (Note that these are apparent 


magnitudes as viewed from the Earth, not from Spitzer). From 
Eqs. (9) and (It) in Ivezic et at. (2001) we can calculate the abso- 
lute, phase corrected magnitudes of Elisa corresponding to these 
apparent values. These magnitudes are H, - 12.29 (prograde) and 
ff r = 12,27 (retrograde). Converting this to V-band absolute magni- 
tudes requires the assumption of a g - r color for Elisa and a trans- 
form equation. The mean g-r color of V-type asteroids is 0,24 
(Moskovitz et al„ 2008). Transform equations can be found in 
Smith et al, (2002) to give phase corrected absolute magnitudes 
of H v = 12.59 (prograde) and Hy* 12.57 (retrograde) for Elisa at 
the time of the Spitzer observations. An average of these absolute 
magnitudes was used in the calculation of Elisa's albedo 
(Section 4). 

The uncertainties associated with these absolute magnitudes 
are dominated by the the assumed g-r color for Elisa. We esti- 
mate that this uncertainty is ±0.04 based on the dispersion of 
g-r colors amongst asteroids in the Vesta family (Moskovitz 
et ai., 2008). There may also be uncertainties introduced by the 
use of a phase correction equation that was derived for S-type 
asteroids (Ivezic et ai., 2001 ). Phase correction equations do not ex- 
ist for V-type asteroids. 

The amplitude of the lightcurve allows an estimation for the ax- 
ial ratio of Elisa. This is achieved by assuming a tri-axial ellipsoid 
and a spin axis aligned perpendicular to the line of sight. Employ- 
ing Eq. (1 ) from Masiero et ai. (2009) and a measured peak-to-peak 
light curve amplitude of 0.35, we estimate Elisa is elongated with 
an axis ratio of at least 1.38. 

3.4. Effects of shape on derived size and temperature 

The thermal models used above to derive the size and temper- 
ature of Elisa assume a spherical shape, but the lightcurve (Fig, 3) 
demonstrates that this asteroid is not a perfect sphere. Departures 
from spherical shape change the surface temperature distribution 
due to different orientations of surface facets relative to the Sun, 
which can affect the radius derived assuming spherical shape. 
We fit the thermal flux spectrum of Elisa with an “ellipsoidal 
STM", assuming axial ratios a:b:c of 1.38:1:1 inferred from the 
lightcurve and the possible rotational phases highlighted in 
Fig. 3. This model makes the same assumptions as the STM, except 
that the temperature solution and integration of flux are for a tri- 
axial ellipsoid rather than a sphere (e.g., Brown, 1985, though our 
implementation using matrix algebra eliminates the need for much 
of the tedious algebra described therein). The resulting best-fit 
effective radius (radius of a sphere with the same projected area 
as the ellipsoid) and sub-solar temperatures were well within the 
uncertainties quoted above for the spherical model (Ar eff < 0.1 km 
and AT mm < 0.5 K). 

We note that the H magnitude at lightcurve mean is slightly 
lower than the H magnitude at the rotational phase of the Spitzer 
observations. Thus, our radius estimate of 5.3 ± 0.4 km at the time 
of the observation corresponds to a radius of 5.2 ± 0.4 km at light- 
curve mean. 


4. Albedo and thermal inertia 


Combining the thermal results with the results of visible pho- 
tometry enables us to solve for the albedo of the asteroid: 


P„ = 7 


Dl 


iD*}(w 04H n 


( 6 ) 


Using the phase corrected absolute magnitude H v from the light- 
curve results in Section 3.3: H v - 12.58 ±0.04, D 0 = 1329.1 ± 
10.4 km (Campins et al„ 1985), and D = 2R= 2(5.3 ± 0.4 km) [Sec- 
tion 3.1, above], we find that the geometric albedo of 956 Elisa is 
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Fig. 4. Smooth-surface thermophysical models of 956 Elisa with equatorial viewing 
geometry. of the best fit to the IRS spectrum of 956 Elisa vs. thermal inertia. The 
best-fit models are 20-30 j m 2 K ' s i in both the prograde and retrograde case. 


Fig. 5. A rough-surface thermophysical model of 956 Elisa, y 1 is minimized for the 
rough-surface model at thermal inertia values of 175 j m 2 K 1 s i {prograde) and 
90 j m 4 K 1 s i (retrograde). 


0.142 ± 0.022. We note that this is significantly lower than the 
albedo of 4 Vesta (0.42 ± 0.05; Tedesco et al, 1992, 2002 f but com- 
parable to that of the near-Earth basaltic Asteroid 3908 Nyx 
(0.16^ 08 /_ 0 ,os; Benner etaL, 2002). 

This albedo estimate may be applied to Eq. (4) by calculating 
the Bond albedo A = qp v (q is the phase integral (Bowell et al., 
1 989) assumed to be 0,39 ± 0.20 since no measurement is avail- 
able), resulting in rj - 1.16 ± 0.05. This value is substantially higher 
than the r/S 0.756 characteristic of large main-belt asteroids 
including 4 Vesta (Spencer, 1990; Lim et al., 2005). 

Since 956 Elisa’s rotation period is roughly three times as long 
as 4 Vesta’s (16.49 h (above) vs. 5.342 h ), if the surface layers of 
the two asteroids were identical in composition and structure 
the effect of thermal inertia would give Elisa the lower p value of 
the two. The thermal inertia of the surface layers of 956 Elisa must 
therefore be greater than that of Vesta's. We therefore conducted 
thermophysical modeling to quantify the constraints on 956 Elisa's 
thermal inertia. 

4. 1. Thermophysical modeling 

Three sets of thermophysical models were calculated in order to 
constrain the thermal inertia of 956 Elisa’s surface layers; a set of 
smooth-surface equatorial cases, a grid of smooth-surface cases 
with varying spin pole geometries, and a rough surface case. 

The first thermophysical case, setting the minimum limit on 
thermal inertia, was calculated for a smooth asteroid surface with 
equatorial viewing geometry (sub-soiar latitude « sub-Spitzer lati- 
tude = 0 15 }, The IRS spectrum of 956 Elisa was then matched to ther- 
mophysical spectra calculated for a grid of thermal inertias from 5 
to 2500 ] nr 2 K" 1 s i for both the prograde and retrograde senses 
of rotation (zero thermal inertia would correspond to the STM 
case); the asteroid’s radius was allowed to float. The results are 
illustrated as Fig. 4. The best-fit thermal inertia solutions for this 
smooth-surface equatorial model are 25-30 J m 2 K 1 s i, corre- 
sponding to « 289.2 K for prograde rotation, and 
T raax = 289.4 K for retrograde. 

Since the spin pole of 956 Elisa is unknown, a similar calculation 
(with a coarser grid of thermal inertias) was then done over a set of 
146 possible spin poles, at steps of 20° in ecliptic longitude and lat- 
itude. No solution produced a thermal inertia below the 


2 A somewhat lower albedo for Vesta (p, a 0.35) may be found using the size 
estimate from Thomas et al. fi997b) but this is still significantly higher than the 
albedo result for 956 Elisa. 



S 10 12 14 16 18 20 22 

Woveiertgth (prr.) 

fig. 6. Spirzer IRS spectrum of 956 Elisa in emissivity. The 13.5-14.15 pm region in 
the IRS, indicated by the grey box, is untrustworthy due to a sporadically appearing 
instrumental artifact [http://ssc.spitzer.caitech.edu/irs/features.htmll. Upper trace: 
SI.1/LL2 scaling of 0,84%, the result if the splice is performed with the unreliable 
"teardrop” region included. Lower trace: SL1/LL2 scaling of 2.57%, the result of 
excluding the "teardrop” region. 


20 J m 2 1C 1 s _ i established in the equatorial model, and nearly 
all fell between 20 and 1 00 j nr 2 K ! s"i The exception is the case 
of poie-on viewing geometry, in which a very high thermal inertia 
(>2500 j nr 2 K _l s‘i a value consistent with solid rock without a 
regolith cover) is possible, but such an observation is unlikely gi- 
ven the 0.4 mag. amplitude of the lightcurve. The distribution 
had r max between 284.3 and 292.2 K, consistent with the 
282.3 ±2.8 K derived from the simpler STM fitting. 

Finally, a rough-surface thermophysical model with equatorial 
viewing geometry was calculated. To present the endmember case, 
100% coverage by hemispherical craters was assumed, correspond- 
ing to RMS surface slope of 66°. Results are shown in Fig. 5. As ex- 
pected, the best-fit thermal inertias are higher than in the smooth- 
surface case: 175 for prograde rotation, 90 J nr 2 K _1 si for 
retrograde. 

In summary, the IRS spectrum of 956 Elisa sets a firm lower lim- 
it of 20 j rrr 2 KC 1 s'i for its surface layer. The most likely thermal 
inertia range is 30-150] nr 2 K' 1 s~i, consistent with unconsoli- 
dated sand or a mixture of rock and fine regolith. Knowledge of Eli- 
sa’s spin pole and sense of rotation would significantly improve 
these estimates. 




6 


LF. Urn et al./ Icarus xxx (2011 ) xxx-xxx 


Table 1 

Sample numbers and size fractions of HED meteorite samples discussed in this work. 


Diogenites 

Ellemeet (coarse) 
Ellemeet (fine) 


0-25 pm 

RELAB MP-D2M-112 
RELAB MP-D2M-113 

LAS* 91 900. 29 


0- 

RELAB TB-RPB-018 

Johnstown 


1000 pm 
0-25 pm 

RELAB MB-TXH-Q9S-A 

Johnstown (Saiisb 

‘ury et al., 1991) 

0-75 pm 

Johnstof.txt 

AlH 77256.96 (Salisbury et a!.. 

0-75 pm 

A77256f.txt 

1991} 

Eucrites 

ALHA 76005 

Basaltic eucrite 

0- 

RE LAB TS-RPB-023 

EET 90020 

Basaltic eucrite 

1000 pm 
0- 

RELAB TB-RPB-020 

PCA 82501 

Basaltic eucrite 

1000 pm 
0-125 pm 

RELAB MP-TXH-124 

ALH 85001 

Cumulate 

0- 

RELAB TB-RP8-01 5 

Yamato 

eucrite 

Cumulate 

looopm 
0-75 pm 

RELAB MT-CMP-061 

980318 

Yamato 

eucrite 
(whole rock) 
Plag, sep. 

0-75 pm 

RELAB MT-CMP-062 

980318 

Yamato 

Px. sep. 

75- 

RE1A8 MT-CMP-063- 

980318 


150 pm 

064 


Table 2 

Synthetic and natural pyroxene samples in Fig. II. 


Synthetic pyroxenes {RELAB; Kiima et al., 2005: Turnock 
etai, 1973) 

RELAB sample number 

Pyroxene composition 

DL-CMP-001 


Etlioo 

Dl-CMP-002 


EltgoFSzo 

DL-CMP-003 


£n?5Fs 5 o 

DL-CMP-004 


EnsoFssc 

BL-CMP-006 


E 1125 FS 73 

DL-CMP-009 


Wo;oEn 4 s Fs4s 

Dt-GMP-010 


Wo te En ra Fs J7 

DL-CMP-012 


Wo s En« s Fs,7.5 

DL-CMP-0I3 


Wo 8 Efts 6 Fs4ii 

DL-CMP-016 


Wo5En 7 S Fs ts 

DL-CMP4M8-A 


Wa,<£n7 2 Fs,s 

DL-CMP-019-A 


W02c£riggFS24 

DL-CMP-021-A 


En 8 fs 52 

DL-CMP-022-A 


Eb?sFs 2 5 

Di,-CMP~025-A 


En 35 F% 5 

DL-CMP-026-A 


En? 0 Fs3 0 

DL-CMP-038-A 


WoaEnsj 

DL-CMP-045-A 


Wd 4 En«Fs48 

ASU pyroxenes (Hamilton,. 2000) 


Sample number 

ASU spectrum 

Pyroxene composition 

BUR-1920 

439 

Woi, 5 £n 77 Fs 2 j 

NMNH- 166555 

541 

Woo. 5 Enr?Fs52 

NMNH-C2368 

549 

Wo4.2En 72 Fs 2 4 

NMNH-1 20414-1 

552 

Wo,, 5 En 7 sFS 22 

NMHH-93527 

558 

Wo 3 En 7 7 Fs 2 j. 


5. Spectral features in the mid-iR 


Wavenumber (cm' 1 ) 

1000.0 666.7 - 500,0 400.0 



Wavelength (>irrs) 


Fig. 7. Spitzer IRS spectrum of 956 Elisa in emissivity, with laboratory spectra 
(RELAB) of diogenites Eliemeet and LAP 91900, cumulate eucrite Y-98G318, and 
basaltic eucrites ALHA 76005 and EET 90020. As in Fig. 6, the grey box represents 
the unreliable 13.5-14.15 pm "teardrop" region of the Sll order of the IRS. Vertical 
dashed lines mark the locations of prominent spectral features discussed in the text, 
such as the “CM" (1 1.4 pm) and "CA7" (19.5 pm) of Hamilton (2000). 


956 Eilsa and Synthetic Pyroxenes 



Woveiength (pm) 

Fig. S. Spectra of 956 Elisa and selected synthetic pyroxene samples plotted in 
Fig. 11. The 13.5-1415 pm region in the IRS. indicated by the grey box, is 
untrustworthy due to an instrumental artifact {Section 3.2) and the transparency 
region (12-13.5 pm; Section 5.1.1) has here been omitted from the laboratory 
spectra for clarity. Vertical dashed lines represent the wavelengths of various 
spectral features. For the 10 pm maximum and the 1 1 .4 pm minimum ("CM”), the 
widths of the lines indicate the formal uncertainties in the central positions from 
Gaussian fits to the spectrum of 956 Elisa. Sample numbers for the pyroxenes are 
given in Table 2. 


In the mid-iR, dividing the spectrum by the thermal continuum 
model yields the spectral emissivity variation of the asteroid 
(Fig. 6) which can then be used to constrain the surface mineralogy 
through comparison with spectra of laboratory samples, including 
HED meteorites, such as are available in the RELAB database 
(http://www.planetary.brown.edu/reiab/) and in Salisbury et al, 
(1991). A list of samples is given in Table 1. 


Comparison with the spectra of HED meteorites (Fig. 7) reveals 
that the closest spectral matches to 956 Elisa are found among the 
diogenites. Like the laboratory spectra of diogenites, the spectrum 
of 956 Elisa exhibits two narrow maxima (10 pm, 11 pm) within 
the broad 9-14 pm silicate band, two maxima at 15.2 pm and 
15,9 pm divided by a 15,55 pm minimum, a minimum in the 
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Fig. 9. Laboratory spectra of coarse and finely-partlculate diogenite sampies. See 
Table 1 for detailed sample identification. Top: RELAB spectra of coarse-grained 
samples from the diogenites Ellemeet (magenta) and LAP 91900 (black). The latter 
is labeled as a 0 to 1000 pm size fraction. Bottom: Spectra Of finely particulate 
samples from three diogenites measured by RELA8 and Salisbury et al. (1991): 
EJletneet (purple; Salisbury et ai„ 1991 ), Johnstown (independent measurements by 
RELAB and Salisbury et al. (1991 )). and ALH 77256 (Salisbury et a!.. 1991). RELAB 
sample numbers MP-D2M-1I3 and MB-TXB-035-A. The two RELAB samples are 
labeled as 0-25 pm size fractions; the two from Salisbury etal. (1991 }are described 
as 0-75 pm particle size. Spectral contrast in the Salisbury et at (1991) samples, 
taken on a biconical reflectance apparatus, has been scaled up by a factor of two to 
enable plotting on the same axes as the RELA8 spectra. RELAB spectra have been 
offset for clarity. Note the reduced spectral contrast and the 12-13-5 pm transpar- 
ency features appearing only in the spectra of the fine particle size samples. 
However, the wavelength positions of the reststrahlen features remain unchanged. 
(For interpretation of the references to colour in this figure legend, the reader is 
referred ro the web version of this article.) 


19-20 (im region, and a Christiansen feature (emissivity maxi- 
mum) just shortward of 9 pm. We discuss these features in more 
detail below, beginning with the 9-14 pm reststrahlen features. 


5. J. Spectral features in the 9- 14 pm reststrahlen band 

Spectral structure within the 9-14 pm reststrahlen band is a 
dose match to that found in laboratory spectra of diogenites. In 
particular, the maxima at 10 and 11 pm and the minimum at 
1 1 .4 pm are narrow, distinct structures that closely match the 
diogenite spectra in wavelength. 

The diogenites are dominated (>84% by volume) by orthopyrox- 
ene of a single very uniform major-element composition 
(WoaMEnwaFSjfci: Mittlefehidt et al., 1998; although a few outli- 
ers extend the Mg range in known diogenitic orthopyroxenes to 
En 8 3 . 8 and En 66 4 (Mittlefehldt et aL, 2009)). We examined the 
spectra of pyroxenes of a range of compositions in order to deter- 
mine how sensitive the spectral features common to diogenites 
and to 956 Elisa are to pyroxene composition. 

The 10- and 11 -pm reststrahlen structures vary in shape and 
wavelength with pyroxene composition (Hamilton, 2000). The 
10 pm emissivity maximum is bounded by the minima designated 
“critical absorptions" CA1 and CA2 by Hamilton (2000). Although 
Hamilton (2000) established the systematic variation of the wave- 
lengths of CA1 and CA2 minima themselves with Fe content, we 
found it impractical to rely on our ability to determine the wave- 
lengths of these minima in the spectrum of 956 Elisa. The short- 
wavelength side of CA1 is governed by the Christiansen feature, 
which may be shifted in wavelength by non-pyroxene minerals 



Fig. 10. Gaussian fit to the 10-pm structure in the emittance spectrum of 956 Elisa. 
The model (magenta trace) is a Gaussian with a linear background; the data are 
represented by the black trace with red error bars. The best-fit Gaussian is centered 
at 9.91 ± 0.05 pm. The same mode! applied to six available diogenite spectra 
resulted in an average Gaussian center of 9.85 1 0.01 pm. (For interpretation of the 
references to colour in this figure legend, the reader is referred to the web version of 
this article.) 


in the rock (for example, olivine) or by pressure and temperature 
effects (Logan et al., 1973). CA2 presented a different problem; at 
the S/N level of the Spitzer data, CA2 and CA3 are not readily dis- 
tinguishable. For these reasons, we instead examined the wave- 
length of the 10 pm emissivity maximum, which also varies 
systematically with the Fe content of low-Ca pyroxenes and was 
readily identifiable in the Elisa data. 

An upper limit on the Mg content of Elisa's pyroxene may be set 
by the absence of the a local minimum at 9.9 pm that divides the 
10 pm structure in magnesian orthopyroxenes En S 3 and above 
(Hamilton, 2000). This also eliminates less magnesian pyroxene 
compositions with Wo > 5% (Fig. 8). The six available mid-lR labora- 
tory spectra of diogenites (Fig. 9) also lack this 9.9 pm minimum; 
the wavelength of the single local maximum in this region falls 
within a narrow range between 9.856 and 9.893 pm irrespective 
of particle size. 

In order to determine the central wavelength and uncertainty of 
this 10 pm structure in the IRS spectrum of 956 Elisa, we applied a 
simple fit consisting of a Gaussian peak on a linear background to 
both the 956 Elisa spectra and the available diogenite laboratory 
spectra (ALH A77256, LAP 91900. johnstown, and coarse and fine 
samples from Ellemeet) from 9.4 to 10.3 pm. The best-fit Gausstans 
for the meteorite spectra are centered around 9.85 ±0.01 pm; the 
best fit to the spectrum ofElisa is 9.91 ±0.05 pm(Fig. 10). The width 
of the best Gaussian fit to the Elisa data is also consistent with the fits 
to the lab data: 0,16 ±0,07 pm for the asteroid vs. 0.20±0.02 pm 
among the laboratory spectra. (The uncertainties represent the for- 
mal uncertainties from the Gaussian fitting.) A similar Gaussian 
model was used to determine the position and uncertainty of the 
11 .4 pm minimum, the ''CA4'' of Hamilton (2000). 

Fig. 11 shows the resulting positions for the 10 pm emissivity 
maxima vs. 1 1 .4 pm minima in the spectra of 956 Elisa and various 
HED meteorites. The ellipses surrounding the points representing 
the diogenites and the basaltic and cumulate eucrites are 2cr from 
their respective means. For comparison, on the same figure, we 
have plotted the positions of these two features in the RELAB spec- 
tra of a collection of synthetic pyroxenes representing a range of 
compositions (Kiima et a!., 2005; Turnock et ai., 1973) and ASU 
emissivity spectra of natural pyroxenes of near-diogenitic compo- 
sition (Hamilton, 2000; Christensen et al., 2000). Diogenites and 
pyroxenes with diogenite-like compositions fall within 2<x of 956 
Elisa, whereas more ferroan, more magnesian, and more calcic 
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Fig. 11. Wavelengths of the 1 1 A fim minimum ("CA^' of Hamilton (2000)} vs. the 10 |xm maximum in the emissjvity spectra of 956 Elisa, and natural pyroxenes of diogenite- 
like composition (purple) from the ASU collection described by Hamilton (2000). Also plotted are the same features as they appear in the inverted reflectance spectra of the 
synthetic pyroxenes (red) of Lmsley and Turnock (RELAB; Klima et ah, 2005, 1973' and of HED meteorites (shown as diamonds; RELAB and Salisbury et ah ( 1991 )). Diogenites 
and pyroxenes with diogenite-hke compositions plot within 2<? of 956 Elisa, whereas more ferroan, more magnesian, and more calcic pyroxenes plot outside this range. See 
Fig. 12. The plotted uncertainties in the location of Elisa’s features are the \u uncertainties emerging from a simple Gaussian fit to each feature (Fig. 10). (For interpretation of 
the references to colour in this figure legend, the reader is referred to the web version of this article.) 


@ Synthetic pyroxenes within 2 “Sigma of 956 E&sa (Erg. \t) 

0 Synthetic pyroxenes more than 2-sigma away from 956 Elisa 
0 Natural pyroxenes .shown in Fig. IX (ASU library) 

# Diogenites with mkf*?R Sab spectra (RElAS or Salisbury) 



F%» 12. Pyroxene quadrilateral illustrating the compositions of diogenites and 
pyroxene samples plotted in fig, 1 1 . 

pyroxenes plot outside this range. We find that the wavelength of 
the CA4 minimum is negatively correlated with Fe content in the 
synthetic pyroxene spectra from Enioo to En, 7 . The compositions 
of the pyroxenes represented on this figure are also plotted on 
the pyroxene quadrilateral in Fig. 12 and their spectra in Fig. 8. 

Pyroxenes in the eucrites are significantly more ferroan than 
diogenitic pyroxenes: consequently, the cumulate eucrites plot be- 
low and to the right of 956 Elisa and the diogenite range on Fig. 11. 
Among the eucrites, basaltic eucrite pyroxenes are more ferroan 
than those of the cumulate eucrites. In Fig. 1 1 , the basaltic eucrites 
consistently appear to the right (longer 10 pm peak) of the cumu- 


late eucrites, but exhibit a wide range of CA4 wavelengths due to 
the broader widths of the CA4 minima in many of the basaltic eu- 
crite spectra. 

Unlike diogenites, eucrites contain substantial amounts of pia- 
gioclase feldspar in addition to their pyroxenes. We were able to 
evaluate the effect of plagioclase on the 1 0 pm and 1 1 .4 pm pyrox- 
ene features in the spectrum of the cumulate eucrite Y-980318, for 
which spectra of both the whole rock and mineral separates 
(plagioclase and two pyroxenes) have been made available in the 
RELAB database (Pieters et al„ 2006). The dominant pyroxene in 
Y-9S0318 is WoioEnsj.s (Buchanan and Tomiyama, 2003), consistent 
with the eucrite’s position in Fig. IT near the synthetic En 50 Fs 5 o 
sample. Notably, the low-Ca pyroxene separates and the whole- 
rock mid-IR spectra have identical wavelength positions for both 
the 10 pm maximum and the 11. 4 pm minimum, suggesting that 
these two spectral features are sensitive to pyroxene composition 
at the plagioclase abundances found in cumulate eucrites. 

Thus, we find that the wavelength positions of the 10 pm max- 
imum and the CA4 1 1.4 pm minimum are diagnostic of pyroxene 
composition, and that in this wavelength region the spectrum of 
956 Elisa is more consistent with a diogenitic (Wo 2 ±T En 7 4 + 2 Fs 24±1 ) 
pyroxene composition than with other HED pyroxenes. 

5.1.1. Transparency feature 

We have not considered the 12-14 pm region in our composi- 
tional mixing models due to the strong effects of particle size in 
this region. Laboratory spectra of fine-grained silicate samples 
exhibit a minimum in emittance, known as the "transparency 
feature," around 12 pm which is absent from the spectra of 
coarse-grained samples (see for example Fig. 9). This feature is also 
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Wavelength (/am) 


Fig. 13. Linear mixing models of the spectrum of 956 Elisa. Vertical dashed bars 
mark the locations of the 10 pm emissivity maximum, the 11 pm maximum, and 
CA4 in the spectrum of Elisa. Top: The first fit considered only the reststrahlen 
features (from 9.0 to 11.7 pm) and resulted in a mixture of approximately 85* 
Ellemeet diogenite and 15* olivine (excluding the spectrally neutral component; 
see Table 3). Bottom*. Fitting region from 7.2 to 11.7 pm. including the Christiansen 
feature. The best fit with the Christiansen feature included 40% Ellemeet diogenite, 
35% cumulate eucrite, 12% olivine, and 12% basaltic eucrite. again excluding a 
neutral component. 

absent from the [RS spectrum of 956 Elisa, which may suggest a 
relatively coarse particle size. 

However, caution is required when interpreting the transpar- 
ency feature on airless bodies as the feature may behave differently 
under vacuum than when measured at atmospheric pressure in the 
laboratory. Henderson et al. (1996) presented emissivity spectra of 
50 pm quartz particles measured in an environment chamber at 
different pressures. A transparency feature between 10 and 
12 pm that is present under higher pressure is absent in the spec- 
trum measured under near-vacuum £1G~ 3 Torr}, Unfortunately, in 
the pioneering study of the effects of pressure on emission spectra, 
Logan et al. (1973) did not measure the transparency region of 
their granite samples. Model calculations by Henderson and Jako- 
sky (1997) for both quartz and basalt reproduce their earlier mea- 
surements, showing a decrease in contrast of the transparency 
feature with decreasing pressure. However, their work indicates 
that the effects of particle size on the transparency feature are 
far more dramatic than the effects of vacuum. The absence of a 
transparency feature in the mid-IR spectrum of Elisa is therefore 
likely to indicate a fairly coarse-grained regolith, but this conclu- 
sion cannot be confirmed without further laboratory data. 

5. 1.2. Linear mixing models 

In order to investigate the presence and abundance of addi- 
tional mineralogical components such as piagioclase, olivine, and 


Wavenumber (cm ') 



1 5 *20 

Wavelength (pm) 


Fig. 14. Linear mixing models from Fig. 13 and endmembers with the spectrum of 
956 Elisa. Spectra of other endmembers, such as Ellemeet. LAP 91900, Y-980318, 
and ALH A76005 may be found in Figs. 7 and 9. 


Table 3 

Results of linear mixing models. 


Component 

Type 

Wavelength range 
9-11.7 pm 7.2-11.7 pm 

Ellemeet (coarse) 

Diogenite 

0.22 

0.20 

Ellemeet (fine) 

Diogenite 


- 

LAP 91900 

Diogenite 

- 

- 

Chassigny olivine 

Fo 69 

0.04 

0.06 

Y-98Q3JQ whole-rock 

Cumulate eucrite 

- 

0.17 

Y-980310 piagioclase 

An 88-91 

- 

- 

Y-980318 Fe-Px 

Wo 2 Fs 48-49 

_ 

- 

Y-980318 Ca-Px 

Wo 45-46 F5 IS- 19 

- 

- 

EET 90020 

Basaltic eucrite 


0.06 

Al H A76005. PC A 82501 
Spectrally neutral 

Basaltic eucrites 

0.74 

0.51 


additional HED materials, we applied a linear mixing model to 
the 7-12 pm region of the Spitzer spectrum of 956 Elisa. As end- 
members, we used three diogenite spectra, three basaltic eucrite 
spectra, the whole-rock spectrum of the cumulate eucrite Y-980318, 
piagioclase (An 88 _.gi) and pyroxene separates from Y-980318, and 
a spectrum of a Chassigny olivine sample Fo 69 , representing 
the olivine composition found in diogenites (Mittlefehldt. 1994) 
and in the HED-related dunite MIL 03443 (Mittlefehldt. 2008), 
Foes-Foys- Although RELA8 spectra, which are measured in 
reflectance at terrestrial atmospheric pressure, are not ideal for 
permitting rigorous quantitative deconvolution of thermal emis- 
sion spectra, these currently represent the best available spectra 
of samples of the appropriate particle sizes for the regolith of 
956 Elisa. Thus, deconvolution against these spectra provides a 
useful qualitative check on which components are most likely to 
be present and a rough estimate of their relative abundances. 

Models of the Christiansen feature (<9pm emissivity maxi- 
mum) using spectral endmembers measured under atmospheric 
pressure must be interpreted with caution. The few mid-IR spec- 
troscopy studies conducted under vacuum (Logan et al., 1973; 
Henderson et al., 1996) have found that the Christiansen feature 
migrates in wavelength at very low pressures (see also Henderson 
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Fig. IS. Spectrum of 956 Elisa from 12.S to 16.S pm. As in Fig. 7, the 13.5-14.15 pm 
region in the spectrum of Elisa is marked by a grey box to indicate the unreliability 
of the IRS in this range. Also shown are the RELAB spectra of two diogenites 
(Eliemeet and LAP 91900), two cumulate eucrites (Y-9803I8 and ALH 85001), and 
two basaltic eucrites (ALH A76O05 and EET 90020). The wavelengths of the 1 5.2 pm 
maximum and 15.55 pm minimum (marked by vertical dashed lines) are a dose 
match in the spectra of Elisa and the two diogenites. 


and Jakosky (1997) in which basalt features were modeled but not 
measured). However, since the experiments both involved very si- 
licic samples (granite and quartz powders), we cannot confidently 
quantify the direction or degree of wavelength change of the Chris- 
tiansen feature of any mafic mineral under vacuum. 

We therefore considered two wavelength ranges: 9.0-1 1.7 pm, 
which isolates the reststrahlen features from the Christiansen and 
transparency regions; and 7,2-11.7 pm, which includes the Chris- 
tiansen feature. Figs. 13 and 14 illustrate the results of the linear 
mixing models; coefficients are given in Table 3. 

When only the reststrahlen features are considered, the best fit 
is given by a combination of the Eliemeet diogenite spectrum at 
coarse (0-1 mm) particle size and olivine, with the diogenite 
approximately five times as abundant as the olivine. No other com- 
ponent improved the fit. This suggests an olivine diogenite compo- 
sition for 956 Elisa. 

However, when the Christiansen feature is included, the best fit 
includes two eucritit components; 40% Eliemeet diogenite, 35% 
cumulate eucrite, 12% olivine, and 12% basaltic eucrite, excluding 
a spectrally neutral component. This model achieves a good match 
to the Christiansen feature at the expense of the CA4 region 
(1 1.4 pm); every component in this mixture apart from the dioge- 
nite has the effect of shifting the model’s CA4 minimum further 
shortward of the CA4 location in the asteroid spectrum. This prob- 
lem is also present to a lesser degree in the reststrahlen-only fit, 
since the olivine in that case influences CA4 in the same way. Thus, 
although these models provide the best overall fits to their respec- 
tive fitting regions, neither matches the CA4 feature as well as a 
single diogenite spectrum. 


Fig. 16. Spectra of plagiociases and pyroxenes from 12.5 to 16.8 pm. Dash-dotted 
portions of the synthetic pyroxene spectra indicate the transparency feature, found 
only 1 in spectra of finely particulate samples (hence absent here from the A5U 
spectra of similar mineralogy). The wavelengths of the 1 5,2 pm maximum and 
1 5.55 pm minimum as found in the spectra of 956 Elisa and diogenites (Fig. 1 5) are 
marked by vertical dashed lines. The 15.55 pm minimum is a weak shoulder in the 
synthetic orthopyroxene spectra (six lower traces). However, the wavelength of the 
15.2 pm maximum is evidently correlated with Mg/(Mg + Fe) composition. Spectra 
have been offset for clarity. 



Tio T95 ZOO 2.05 

Bil Center 


Fig. 17. NiR band centers of 956 Elisa compared with those of 4 Vesta and the HED 
meteorites. The errors on the band centers for the asteroids are comparable to the 
size of the plotted symbols. By direct comparison to the meteorites, the band 
centers of Elisa are most similar to the howardites. 


Other components in the library (Table 3), including the LAP 
91900 diogenite, were excluded from both fits. Eliemeet may pro- 
vide a better match than LAP 91900 because of its greater otivine 
content; Eliemeet has 4.7% modal olivine, whereas LAP 91900 has 
under 1% (Bowman et al„ 1997). Combining EHemeet’s olivine with 
the olivine from the Chassigny component in the linear mixture 
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Fig. 18. Modified Gaussian Model (MGM) fits to She NiR spectrum of 956 Elisa, 
MGM band parameters are consistent with a composition for Elisa that is 
approximately 75* magnesian orthopyroxene (e.g. diogenite-like) and 25% high- 
Ca pyroxene. The average HCP/(HCP+LCP) value from the MGM fits is 26% with 
variation of order ±5%. 


raises the olivine proportions in the mixtures to 18% and 1 4% for the 
reststrahlen-only and Christiansen-feature models respectively. 

Apart from the uncertainty in the behavior of the Christiansen 
feature itself under vacuum, another possible explanation for the 
compositional differences between the two regions is particle size: 
unlike that of the Christiansen feature, the spectral contrast of 
reststrahlen features is known to be reduced at fine particle sizes 
(e.g, Logan et ai., 1973), Thus, the Christiansen feature is more 
likely to reflect the bulk composition of the material, whereas 
the reststrahlen features will strongly overrepresent the coarsest 
particulate material in the field of view. 

5.2. Spectral features, 14-17 pm 

Hamilton (2000) notes that features in the 800-600 cm" 1 re- 
gion (12.5-16.7 pm) are produced by symmetric Si-O-Si stretch- 
ing vibrations and weaken in band strength with increasing Fe, 
but does not evaluate their diagnostic potential because of their 
low spectral contrast in the ASU library pyroxene samples. 

However, the spectrum of 956 Elisa exhibits maxima at 15.2 
and 15.9 pm. separated by a 15.55 pm emissivity minimum, that 
are closely matched in wavelength by the RELAB spectra of dioge- 
nites (Fig. 15). In cumulate eucrite spectra in this region a similar 
feature occurs at a slightly longer wavelength; in basaltic eucrite 
spectra the minimum is significantly weaker and occurs at a longer 
wavelength still. 

In the RELAB spectra of synthetic iow-Ca pyroxenes (Fig. 16), 
the wavelength of the 15.2 pm maximum is evidently correlated 
with En content. The wavelengths of this feature in the spectrum 
of the En g o and Wo.,En-?6 samples are the closest matches to its 
wavelength in the diogenite spectra (and in the spectrum of 956 
Elisa), a result consistent with the known diogenite pyroxene 
compositions. 

We note that the two diogenite spectra in Fig. 15 also display a 
narrow minimum at 14.6 pm which is not evident in the spectrum 
of 956 Elisa; however, it is also absent from spectra of several of 
the finely-particulate diogenite samples (Fig. 9). 

5.3. 16-25 pm reststrahlen features 

Four “critical absorptions” were identified in orthopyroxenes by 
Hamilton (2000) in the 16-25 pm reststrahlen region at 17.6. 18.4, 
19.5, and 22.1 pm, with the observation that none of them deviates 


Table 4 

MGM band parameters. 



Center (pm) 

Width (pm) 

Strength (log reflectance) 

4-JuIy-2008 



Band 1 

0.277 

0.223 

-0.40 

Band 2 

0.376 

0.264 

-0.57 

Band 3 

0:629 

0.219 

-0.28 

Band 4 

0.902 

0.216 

-0.75 

Band 5 

0.995 

0.189 

-0.23 

Sand 6 

1.173 

0.277 

-0.16 

Sand 7 

1.897 

0.602 

-0.62 

Sand S 

2.268 

0.587 

-0.16 

Sand 9 

2.600 

0.505 

-8E-4 

Continuum intercept * 0.945, slope* -0.021 prrT : 
Residual « 1.20% 

S -October-2008 


Band 1 

0.280 

0.289 

-0.74 

Band 2 

0.374 

0.278 

-0.52 

Band 3 

0.626 

0.237 

-0.36 

Sand 4 

0.902 

0.219 

-0.78 

Band 5 

0.990 

0.190 

-0,21 

Band 6 

1.172 

0.282 

-0.14 

Band 7 

1.906 

0.573 

-0.67 

Band 8 

2.262 

0.590 

-0.12 

Band 9 

2.603 

0.516 

-0.02 

Continuum intercept - 0.882, slope- -0.013 pm' 1 
Residual -1.11% 



substantially in wavelength as a function of Mg or Fe content. Of 
these, only the deepest (CA7, 19.5 pm) is readily identifiable in 
the spectrum of 956 Elisa. CA7 is found at the same wavelength 
in the spectra of diogenites; at a lower spectral contrast and 
slightly longer wavelength in the spectra of cumulate eucrites; 
and not at all in the spectra of basaltic eucrites. 


6. Near-IR spectroscopy 

Near-infrared (NIR, 0.75-2.5 pm) spectra of 956 Elisa were ob- 
tained to complement the mid-IR Spitzer observations. Data were 
obtained on 4-july-2008 and 9-0ctober-2008 at NASA’s Infrared 
Telescope Facility (1RTF) using the SpeX instrument (Rayner 
et al, 2003), Reduction of these data employed the IDL-based Spex- 
tool package (Cushing et al., 2004). Since these data did not display 
significant spectral variability as a function of Elisa’s rotationai 
phase (the longest base line that was tested was —3 h on 9-Octo- 
ber-2008), we combined individual exposures from each night into 
two final spectra. The maximum time for individual exposures was 
120 s and the net exposure time for 4-july-2008 was 16 min and 
for 9-October-20Q8 was 88 min. Visible-wavelength data from 
Lazzaro et al. (2004) were appended to our N!R data to extend 
the spectral range down to 0.5 pm. 

A simple band analysis in the style of Cloutis et al. ( 1 986) was 
performed for each of the NIR spectra. Specific definitions of the 
band parameters measured here can be found in Moskovitz et al. 
(2010), The temperature-corrected (Burbine et al., 2009) 1 and 
2 pm (BI and BH respectively) band centers were measured at 
0.932 and 1.962 pm for the data from 4-ju!y-2008 and at 0.931 
and 1 .938 pm for the data from 9-October-2008. 

The temperature corrections were calculated using the mea- 
sured albedo of 0.142 (Section 4) and were found to be 0.002 pm 
for BI and 0.01 7 pm for Bli, These corrections are the same for both 
nights because the heliocentric distance of Elisa happened to be 
the same. The error bars on these band centers are no greater than 
0.002 pm for BI and 0.01 for BH. These band centers are plotted in 
Fig. 17 along with those of asteroid 4 Vesta and the regions occu- 
pied by each of the HED meteorite subgroups (Moskovitz et ai., 
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201 0), This figure shows that by spectral analogy the band centers 
of Elisa are most consistent with a howardite-like composition. 
However, it is important to note that these band centers represent 
a simplified parameterization of the 1 and 2 pm absorption fea- 
tures. in reality these features are produced by several overlapping 
bands, each with distinct central wavelengths. As a result the mea- 
sured band centers merely provide approximate indications of bulk 
composition. 

The measured band centers can be used with the equations of 
Caffey et al. (2002) and Burbine et al. (2009) to compute bulk sin- 
gle-pyroxene Fs and Wo molar contents. As the difference in mea- 
sured band centers from the two nights of observing are not 
significant in these calculations, average values are presented. 
The equations of Gaffey et al. (2002) predict Fs = 40 + 5 and Wo 
= 9 + 4. The equations of Burbine et a!. (2009) predict values of 
Fs = 38 + 3 and Wo = 8 ± 1. These Fs-Wo numbers suggest a pyrox- 
ene composition for Elisa that is consistent with polymict eucrites 
and howardites (Burbine et al.. 2001, 2009). 

These results merit further investigation, particularly in light of 
our mid-lR analysis (Section 5) which suggests a composition dom- 
inated by diogenite-like, low-Ca pyroxene (ICP). Thus, we employ 
the modified gaussian model (MGM; Sunshine et a!., 1990) to pro- 
vide a more realistic treatment of Elisa's NIR absorption features. 
Fig. 18 shows the results of this MGM analysis and Table 4 summa- 
rizes the parameters of the fitted bands. Following convention (e.g. 
Sunshine et al., 1990) the data have been renormalized so that the 
peak reflectance around 1 .4 pm is less than unity, hence the nega- 
tive values of ln( reflectance) shown in Fig. 18. In short, we were 
unable to fit Elisa’s NIR spectrum using only LCP: both LCP and 
HCP were required to produce a reasonable fit. 

Nine bands were used to model the reflectance spectra of Elisa: 
three visible wavelength bands centered around 0.3, 0.4 and 
0,6 pm , the primary LCP bands near 0.9 and 1.9 pm , the primary 
HCP bands near 1.0 and 2.2 pm , one band around 1.2 pm and an 
atmospheric water band around 2.6 pm. The modeling of pyrox- 
ene-dominated compositions with these nine bands is consistent 
with previous studies (Sunshine et al., 2004, Mayne et al., submit- 
ted for publication). The residual root-mean-square (RMS) errors 
on these fits were ~1%. 

The MGM bands display similar central wavelengths, widths 
and strengths in the two spectra. The variation that does exist (Ta- 
ble 4) is most likely a consequence of bands shifting to compensate 
for differences in the residual atmospheric features near 1.4 and 
1.9 pm. Comparison of Elisa’s LCP band centers (Bands 4 and 7 in 
Table 4) to those of synthetic pyroxenes (Klima et al., 2007) suggest 
a Mg# certainly greater than 50 and probably closer to ~80. The ra- 
tio of LCP to HCP band strengths in both the 1 and 2 pm bands give 
an indication of the fraction of surface material composed of HCP 
(Sunshine and Pieters, 1993). The average HCP/(HCP+ LCP) value 
from the MGM fits is 26% with variation of order ±5%. 

In total, these MGM band parameters suggest a composition for 
Elisa that is approximately 75% magnesian orthopyroxene (e.g. 
diogenite-like) and 25% high-Ca pyroxene. This is consistent with 
the result of our simple band analysis (Fig. 17). The presence of 
high-Ca pyroxene pushes Elisa's R! and Bil centers towards longer 
wavelengths than would be consistent with a composition of 
exclusively diogenite-like material. Considering that the band 
analysis performed for Fig. 1 7 provides a rough, first-order approx- 
imation of mineralogy, we feel that these two techniques provide 
complementary insight on the composition of Elisa. 


7. Discussion and comparison with near-IR spectra 

Based on the wavelengths of the 10 pm maximum, the 1 1 .4 pm 
minimum, and the 15.2-1 5.9 pm Si-O-Si symmetric stretching 


features, as well as the strength of the 1 9.5 pm minimum, we con- 
clude that the Spitzer IRS spectrum of 956 Elisa is most consistent 
with thermal emission from a material of predominantly oiivine- 
diogenitic composition. Non-cumulate (basaltic) eucrite spectra 
cannot account for the observed features. The general pattern of 
reststrahien features in the cumulate eucrite spectra is similar to 
that found in diogenites and the IRS data, but the CA4 minimum 
(11.4 pm), 15.2-15.9 pm features, and CA7 minimum (19.5 pm) 
are all somewhat poorer matches to 956 Elisa in wavelength than 
are the diogenite equivalents. Comparison with the library suites of 
pyroxenes establishes that these small shifts in wavelength are 
compositionally meaningful. 

High-Ca pyroxene is present at 25% ± 5% in the MGM solution to 
the near-IR spectra. Two-pyroxene solutions to the spectra of Ves- 
toids have been uncommon in the literature. This is a considerably 
lower level than the 59-66% found for the three V-type near-Earth 
asteroids (3908 Nyx, 4055 Magellan, and (6611) 1993 VW) by Ca- 
ms et al (2008), and slightly lower than the 30-40% found for 
41 88 Kitezh by Sunshine et al, (2004). Calcic pyroxene is not evi- 
dent in the mid-IR spectrum of 956 Elisa, but may belong to a com- 
ponent that is underrepresented in the mid-IR spectrum either 
because of its spatial distribution on the asteroid or because of 
its particle size. In addition, high-Ca pyroxene may be present as 
exsolution lamellae: fine-scale exsolution has been found in some 
diogenites, and the natural orthopyroxene samples of Hamilton 
(2000) were sometimes noted to have included exsolution impuri- 
ties that appeared not to have affected the mid-IR spectra. Thus, it 
is possible that exsolved pyroxenes may be less evident in thermal- 
IR spectra than they are in the near-IR. 

Among HED materials, high-Ca pyroxene occurs in eucrites and 
is usually accompanied by plagioclase. Although various eucrites 
and a eucritic plagioclase were included in the spectral library 
for the fitting of the 9-11.7 pm reststrahien features, the best fits 
to this region were plagioclase-free. In contrast, the fits to the 
Christiansen feature consistently produced compositions contain- 
ing substantial amounts of eucritic material. This combination sug- 
gests a particle size effect: at small particle sizes, spectral contrast 
in the reststrahien features is suppressed, but spectral contrast in 
the Christiansen feature is not. 

The hypothesis that eucritic and diogenitic material on 956 Eli- 
sa occur at different particle sizes suffers from the obvious draw- 
back that the two lithologies do not have very different 
mechanical properties. One possible solution is that a recent im- 
pact into a diogenitic terrane on Elisa may have excavated fresh 
orthopyroxenite. 

Thermal infrared observations are also strongly biased toward 
the warmest part of the asteroid in comparison with near-IR obser- 
vations. For example, 76% of the IRS flux would have originated 
within 50° of the sub-solar point, a region representing 37% of 
the visible hemisphere. Nevertheless, the IRS and SpeX observa- 
tions are both consistent with the presence of a substantial amount 
of diogenitic pyroxene observable on 956 Elisa, 

In summary, on the basis of Spitzer IRS data and new ground- 
based observations of 956 Elisa, we determine that the rotation 
period of 956 Elisa is 16.494 ±0.001 h; that its geometric albedo 
is p B = 0.142 ± 0.022; that its radius is 5.3 ± 0.4 km; that its maxi- 
mum surface temperature at the time of observation was 
282.3 ± 2.8 K; and that its thermal spectrum is a good match to 
the standard thermal model with = 1.16 + 0.05. Thermophysical 
modeling of the IRS spectrum of Elisa places a firm lower limit of 
20 j nrr 2 K“ ! s4 on the thermal inertia of the asteroid’s surface 
layer (if the surface is very smooth) but more likely values fall be- 
tween 30 and 150 j nr 2 K' ! s~i depending on the sense of rota- 
tion, which remains unknown. 

Spectral features observed within the 9-12 pm reststrahien 
band, the 15-16,5 pm symmetric Si-O-Si region (Hamilton, 
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2000), and the 16-25 jim reststrahlen region are closely matched 
by features found in laboratory spectra of diogenite meteorites. 
Non-cumulate (basaltic) eucrite spectra cannot account for the fea- 
tures observed in the IRS spectrum of 956 Elisa; linear spectral 
deconvolution against a suite of BED meteorite spectra and related 
samples suggests that they can be excluded as a major component 
of the portion of 956 Elisa represented by the IRS spectrum. Cumu- 
late eucrite material may be present at lower abundances than the 
diogenitic component. However, the diagnostic 9-12 pm reststrah- 
len features are better fit by a mixture of diogenite and 10-15% 
olivine (e.g. olivine diogenite) than by a mixture that includes 
cumulate eucrites. 
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